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The phonon contribution to the shear viscosity 77 in superfluid neutron stars is calculated by 
assuming neutron pairing in a 'So channel. The shear viscosity is obtained by means of variational 
methods for the solution of the Boltzmann equation amended by a collision term which takes into 
account the binary collisions of phonons. Effective field theory techniques are used to extract 
the phonon scattering rates in terms of the equation of state (EoS) of the system. We find that 
77 oc 1/r 5 , the proportionality factor depending on the EoS of the system. Our results indicate 
that the phonon contribution to 77 might have important effects for the different oscillation modes 
of the star. 
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1. Introduction 

The study of superfluidity in compact stars, such as neutron stars, has received a lot of attention 
after Migdal's observation that superfluidity of neutron matter may occur in the core of compact 
stars In fact, this possibility has been taken into account to explain very different neutron star 
phenomena, such as pulsar glitches, star cooling or the dynamics of the neutron star's oscillations. 

At low temperatures neutron matter superfluidity occurs after the appearance of a quantum 
condensate, associated to neutron pairing. The condensate spontaneously breaks the global U(l) 
symmetry associated to baryon number conservation. Thus, the Goldstone's theorem predicts the 
existence of a low energy mode which is essential to explain the property of superfluidity, the 
so-called superfluid phonon. 

In this manuscript we will explore the phonon contribution to the shear viscosity rj of super- 
fluid neutron matter in the core of the stars, assuming that the neutrons pair in a So channel [0]. In 
order to compute the phonon contribution one needs to asses the relevant phonon collisions which 
are responsible for the transport phenomena. The leading phonon interactions can be completely 
determined by the equation of state (EoS) of the system [[|^|]. In this paper we exploit the universal 
character of the effective field theory (EFT) at leading order to present a very general formulation, 
that only depends on the equation-of-state (EoS) of the system. 

2. Phonon interactions in superfluid matter 

The superfluid phonon is the Goldstone mode associated to the spontaneous symmetry break- 
ing of a U(l) symmetry, which corresponds to particle number conservation. EFT techniques can 
be used to write down the effective Lagrangian associated to the superfluid phonon. The effective 
lagrangian is given by an expansion in derivatives of the Goldstone field, the terms of this expansion 
being restricted by symmetiy considerations. The coefficients of the Lagrangian can be computed 
from the microscopic theory, through a standard matching procedure, thus depending on the EoS 
of the system. 

At the lowest order in a derivative expansion, the Lagrangian reads [|[ 

JSf L o = P(X) , 

X = n-d t q>- i-^-, (2.1) 
2m 

where P(ju) and jU are the pressure and chemical potential, respectively, of the superfluid at T = 0. 
The quantity (p is the phonon field and m is the mass of the particles that condense. After a Legendre 
transform, one can get the associated Hamiltonian, which has the same form as the one used by 
Landau to obtain the self-interactions of the phonons of 4 He [Q, ||]. 

After a Taylor expansion of the pressure, and rescaling of the phonon field to have a canoni- 
cally normalized kinetic term, the Lagrangian for the phonon field is given by 

^LO =l((dr$) 2 -v 2 ph (V<P) 2 ) 

-g {(d^) 3 - 3T] g d f <HV0) 2 ) + A ((d f 0) 4 - ifo,! (^0) 2 (V</.) 2 + T] A , 2 (V0) 4 ) + • • • (2.2) 
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The different phonon self-couplings of Eq. ( |Z2| ) can be expressed as different ratios of derivatives 
of the pressure with respect to the chemical potential [§], or in terms of the density, the speed of 
sound at T = 0, and derivatives of the speed of sound with respect to the density. The speed of 
sound at T = is given by 




(2.3) 



where p is the mass density, related to the particle density p as p = mp. 

The dispersion law obtained from this Lagrangian at tree level is exactly E p = c s p. And the 
three and four phonon self-coupling constants can be expressed as 

- 1 (x-7^\, ng - c - 



6-y/mp c s \ c s dp J ' 8 6y/mp g ' 

1 ^l-8^^ + 10^^V-2- p2 ^ 



24 mp c 2 I "c s dp ' '"cj \dp J c s dp 2 J 

n*i= r4t' ^= iT ¥- (2 - 4) 

8 mp A r\ g 

3. Equation of state of superfluid neutron matter 

The speed of sound at T = as well as the different phonon self-couplings are determined by 
means of the EoS for neutron matter in neutron stars. A common benchmark for a nucleonic EoS is 
the one by Akmal, Pandharipande and Ravenhall (APR) [Q] in /3 -stable nuclear matter. Heiselberg 
and Hjorth- Jensen parametrized the APR EoS of nuclear matter in a simple form [0], which will 
subsequently be used in this paper. The effect of neutron pairing is not considered because it is not 
expected to have a big impact in the EoS. 

Using this EoS, one can obtain the speed of sound and all three and four phonon self-couplings. 
The speed of sound can be computed as 



1 dP N 

(3.1) 



m dp„ ' 



where to compute the mass density we only take into account the nucleonic part to the pressure (Pat), 
asm> m e , with m and m e being the mass of the nucleons and electrons, respectively. Further, in 
j3 -equilibrated matter, p « p n , where p„ is the neutron density. For the phonon self-couplings, we 



will substitute p — > p n in Eq. (pL4|). 

The ratio of the speed of sound c s with respect to the speed of light c is shown in Fig. [I] for 
both a free gas of neutrons and for j8 -stable nuclear matter as a function of the density. We observe 
that the assumption that the APR model for /3 -stable nuclear matter is non-relativistic breaks down 
at densities of the order of 1.5-2 po- For those densities, relativistic effects appear as the APR EoS 
includes not only two nucleon but also three nucleon interactions. This seems to suggest that at 
very high densities we should employ the relativistic version of the phonon effective field theory 
[Q]. This would imply, though, very minor changes in the computation of r\ (see Ref. [Q]). 
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4. Shear viscosity due to the binary collisions of superfluid phonons 

The shear viscosity T7 is as a dissipative term in the energy-momentum tensor Ty. For small 
deviations from equilibrium one finds that 

STij = —TjVij = —Vj (diVj + djVi - |fyV- VJ , (4. 1) 

where V is the three fluid velocity of the normal component of the system. 

The superfluid phonon contribution to the energy-momentum of the system is given by 



2 f d3 P PiPj 



where / is the phonon distribution function. The distribution function obeys the Boltzmann equa- 
tion [§ 

being in the superfluid rest frame, and C[f] is the collision term. For the computation of the 
collision term it is enough to consider binary collisions. This is due to the fact that the next-to- 
leading corrections to the dispersion relation of the phonons in superfluid neutron matter implies 
that the decay of two phonons into two is the first kinematically allowed process [g]. 

In order to compute the shear viscosity we proceed by considering small departures from 
equilibrium to the phonon distribution function and linearizing the corresponding transport equa- 
tion [Q]. It is then necessary to use variational methods in order to solve the transport equation, as 
in Refs. [^, [R], 11]. The final expression for the shear viscosity is given by [Q] 



V= \T5j TfM' (4 " 4) 

where M is the is the 2f>2 scattering matrix for phonons. After a very simple dimensional analysis, 
we observe that the shear viscosity due to the binary collisions of superfluid phonons scales with 
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Figure 2: Left: Shear viscosity multiplied by T 5 (rj T 5 ) as function of density for a free neutron gas and 
j3 -stable nuclear matter. Right: Phonon mean-free path in /3 -stable matter as function of temperature for 
different densities. The straight line indicates the hydrodynamic limit: Ih < 10 Km (taken from Ref. ||^]). 



the temperature as l/T 5 , a universal feature that occurs in other superfluid systems such as He 
or superfluid cold atoms at unitary [[k], 12j |. 



5. Discussion 

We start this section by studying the p dependence of 77 on the left-hand side of Fig. || as 
we exploit that r\ T 5 is a temperature-independent quantity. We consider both the APR EoS that 
describes j3 -stable nuclear matter and the EoS of a free neutron gas. The viscosity increases by 
almost three orders of magnitude when the density changes from po to 2 po for /3 -stable matter, 
the increase being much more moderate for a free neutron gas. This is related to the behavior of c s 
seen in Fig. |l|. We can then conclude from Fig. || that the different choices of EoS have a very clear 
impact in the numerical values of r\ while there is a universal l/T 5 behaviour with temperature. 
Needless to say, one should also study other processes such as the phonons colliding with other 
particles in the star, such as the electrons, in order to assess the dominant processes for dissipation 
in stars. We leave this computation for future studies. 

A relevant discussion also is to know the temperature regime where our results are applicable. 
Hydrodynamics is only valid when the mean free path is smaller that the typical macroscopic length 
of the system, in this case, the radius of the star. We show on the right-hand side of Fig. |2| the mean 
free path of phonons within j3 -stable matter for 0.5po, po and 2p as a function of temperature. We 
also indicate an estimate of the limit of 10 Km, the radius of the star. The mean free path / was 



extracted from the computation of rj in Ref. flllfl 



/ 



'? 



n < p > 

where < p > is the thermal average momentum, and n is the phonon density: 

T 



(5.1) 



< p >=2.7- 



d 3 p T 3 



(27Cf 



2„3 ' 



7I A C 



(5.2) 
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On the right-hand side of Fig. || we observe that for temperatures below T ~ 0.1 MeV, the 
phonon mean free path is bigger than the size of the star. This would seem to indicate that for 
T < 0. 1 MeV it is questionable to have a hydrodynamical description of the phonons of the star. 
We note that the critical temperature for the phase transition to the normal phase is T c ~ 1 MeV, 
although this value, as well as the values we obtained for the phonon mean free path, are dependent 
on the EoS. 
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